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I. INTRODUCTION 



This paper continues our study of nonzero ground state entropy, So({G},q) ^ 0, i.e., ground state de- 
generacy (per site) W({G},q) > 1, where Sq = fcslnVF, in q-state Potts antiferromagnets on various 
lattices and, more generally, families of graphs {G}. There is an interesting connection with graph theory 
here, since the zero-temperature partition function of the above-mentioned g-state Potts antiferromagnet on 
a graph G satisfies Z(G 7 q, T — 0) paf — P{G 1 q), where P(G, q) is the chromatic polynomial expressing the 
number of ways of coloring the vertices of the graph G with q colors such that no two adjacent vertices have 
the same color ||- ||. Thus, 

W([\im G],q) = lim P{G,q) 1/n (1.1) 

n — >oo n — >oo 

where n = v(G) is the number of vertices of G. An example of a real substance exhibiting nonzero ground 
state entropy is ice (t]||. Just as complex analysis provides deeper insights into real analysis in mathematics, 
the generalization from q G Z+, to q e C yields a deeper understanding of the behavior of W({G},q) for 
physical (positive integral) q. In general, W({G},q) is an analytic function in the q plane except along a 
certain continuous locus of points, which we denote B. In the limit as n — > oo, the locus B forms by means 
of a coalescence of a subset of the zeros of P(G, q) (called chromatic zeros of G). Q In a series of papers we 
have calculated and analyzed iy({G}, q) for a variety of families of graphs Jll[]- both for physical values 
of q (via rigorous upper and lower bounds, large-g series calculations, and Monte Carlo measurements) and 
for the generalization to complex values of q. 

A basic question that one can ask about the locus B for (the n — ► oo limit of a) family of graphs {G} 
is whether it extends only a finite distance from the origin of the q plane or, instead, extends an infinite 
distance away from this origin (i.e., passes through the point z = in the z = 1/q plane), and hence is 
noncompact. The importance of this question for the statistical mechanics of q-state Potts antiferromagnets 
is that large-q Taylor series provide a powerful means of obtaining approximate values of the ground state 
degeneracy even for moderate values of q |l7]]- JhJ, p^JI^ ], but these exist if and only if the reduced function 
Q WredX{G}, q) = q~ l W({G}, q) is analytic at 1/q = 0, which, in turn, is true if and only if the nonanalytic 
boundary B does not pass through z — 0. Large-g series calculations of W re d. have been derived for regular 
lattices Jn|, but we noted jl2| that for the r — > oo limit of the graph K2+C r the locus B, which is noncompact 
|jTof , no such large-g series exists. Here, K p denotes the complement of the complete graph K p Qlt is clearly 
important to understand better the differences between the behavior of Potts antiferromagnets on families 
of graphs that yield W re d. ({G}, q) functions analytic at z = and those that do not, i.e. those that yield 
noncompact loci B passing through z = In the present paper, we carry out such a study for families 
of graphs which are homeomorphic expansions of the following type: we start with families of the form 
(K p )t, + G r , where G r is a given graph family with r vertices, b signifies the removal of certain bonds in the 



1 This reduced function is a natural object to define since an obvious upper bound on P(G, q) describing the 
coloring of an n-vertex graph with q colors is P(G,q) < q n , and hence W({G},q) < q; this guarantees that 
lirriq^oo W re d. {{G} , q) is a finite quantity (even if it is sometimes nonanalytic at this point). The large-q series 
take their simplest form as series in the expansion variable y = l/(q — 1) for the equivalent reduced function 
W({G}, q) = q x {l — q 1 )~ A ^ 2 W({G}, q), where A is the coordination number of the lattice. 

2 The complete graph K n on n vertices is defined as the graph where each of these vertices is connected to all of 
the other n — 1 vertices by bonds. The complement K n of K„ is the graph with n vertices and no bonds. 
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K p subgraph, and G + H is the "join" of the graphs G and H We then perform homeomorphic expansion 
of the bonds connecting the K p and G r subgraphs. This family is categorized as being of type HEC . An 
interesting feature of the families analyzed here is that, in contrast to those that we have previously studied, 
they yield loci B (which are boundaries of regions of analyticity of W) that have support for Re(q) < 
and for which the point z = is a multiple point on the algebraic curve B in the technical terminology 
of algebraic geometry (^0|, i.e., a point where several branches of this curve cross each other. For basic 
definitions on homeomorphic classes of graphs, see Ref. ||; some recent work on homeomorphism classes in 
a different direction than ours is Ref. |p2| . 

A general form for the chromatic polynomial of an n- vertex graph G is 

N a 

P(G,q)^c (q)+Y, C ^l) a ^l) t3n (1-2) 

where Cj(q) and a,j(q) are certain functions of q. Here the dj{q) and Cj^o(q) are independent of n, while 
Co(q) may contain rt-dependent terms, such as (—1)™, but does not grow with n like (const.) 11 . A term 
ag(q) is defined as "leading" if it dominates the n — > oo limit of P(G,q); in particular, if N a > 2, then it 
satisfies |a^(g)| > 1 and |o^(g)| > |aj(g)| for j ^ i, so that \W\ = \a^\ tj . The locus B occurs where there is a 
nonanalytic change in W as the leading terms at in eq. (|l.2| ) changes. 

Since for the families of graphs studied here, the boundary B is noncompact in the q plane, it is often 
more convenient to describe the boundary in the complex z or y planes, where, as above, z = 1/q, and 

1 z 

y = — 7 = 1 — (1.3 

q—l 1 — z 

The variable y is commonly used in large-g series expansions. For the families considered here B is actually 
compact in the y and z planes. We define polar coordinates as 

z = (e ie (1.4) 

and 

V = pe lf} (1.5) 

and the function 

Dk{q) = T^s = a " 2 £ 2( - a) " = £ ( - i)s ( k ~s ') ^ s (L6) 

qyq 1 j=0 s=0 V ' / 

where a = q — 1 = 1/y and P(Ck,q) is the chromatic polynomial for the circuit (cyclic) graph CV- with k 
vertices, 

P(C k ,q) =a k + (-l) k a (1.7) 



3 The "join" of two graphs G and H is obtained by adding bonds connecting each of the vertices of G to those of 
H. This was denoted G x H in Ref. 113]; here we use the more common notation for this object in the mathematical 
literature, viz., G + H . 
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We shall also use a standard notation from combinatorics, 




(1.8) 



The organization of the paper is as follows. In section 2 we construct several general multiparameter 
homeomorphic classes of families of graphs that, in a certain limit, yield noncompact boundaries B(q). In 
sections 3-5 we give exact calculations of the respective boundaries B for three such families. These exhibit 
features going beyond those that we found in our earlier work in three main ways: (i) the point z — can 
be a multiple point; (ii) B includes support for Re(q) < 0, and (hi) it is no longer true in general that 
Re(q) = q = for q S B. A general discussion of our findings is given in section 6, and our conclusions 
are presented in the final section. 



We can construct a large variety of families of graphs with noncompact boundaries B(q) of regions where the 
respective W functions are analytic by performing homeomorphic expansions of the basic family constructed 



where, as before, b denotes the removal of b bonds connecting a vertex v in the K p subgraph to the other 
vertices of K p . Since each vertex v € K p has degree A = p — 1, b is bounded above according to b < p — 1. 
We have shown earlier that (i) the locus B for lim r ^ 00 [(_ft'p)b + G r ] is noncompact in the q plane, passing 
through z = 1/q = J51| and (ii) the analogous locus for the r — > oo limit of homeomorphic expansions of 
this family has the same noncompactness property. To construct specific homeomorphic expansions, select 
a vertex v± in the K p subgraph and perform homeomorphic expansion on each of the bonds connecting this 
vertex with the vertices in G r by inserting k\ — 2 additional degree-2 vertices on each of these bonds, where 
ki > 3. Then continue this process with a second vertex V2 € K p , inserting hi — 2 additional degree-2 vertices 
on each of bonds connecting vi with the vertices of G r , and so forth for the other vertices in K p (including 
the vertex v from which the b bonds were removed, if there are any bonds left linking with this vertex). We 
denote the resultant homeomorphic expansion as 



The labelling is chosen so that, counting the two vertices on the original bond connecting to Vj € K p together 
with the kj — 2 inserted vertices, there is a total of kj vertices in all along what was originally this single 
bond. Some illustrative examples of this and other types of families to be discussed are shown in Fig. [l| 



II. HOMEOMORPHIC CLASSES OF HEC TYPE 



and analyzed before [H| , 



(K p ) b + G 



(2.1) 



HEC kl - 2 2[(Kp)b(v) + G r ] 



(2.2) 
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(c) (d) _ _ 

FIG. 1. Illustrations of families of graphs considered in the text: (a) H kr = H ECk 1 -2,k 2 — 2(^2 + K r ), with r = 4 
and k\ = 3, fc 2 = 4, whence k = ki + fc 2 — 1 = 6; (b) &k lt r = HECk 1 -2,o,o(^3 + ^V) with fci = 3 and r = 4; (c) 
Ofc 1>r = HEC kl _2,o(l<2 + TV) with fci = 3 and r = 4; (d) l/ fcj) . = HEC kl -2,k 2 -^2 + T r ) with fci = fc 2 = k = 3 and r = 4. 
Here, K p is the complete graph with p vertices and T r is the tree graph with r vertices. 

Another starting family with a noncompact B in the r-n» limit is 

(-Kp){b};M + G r (2.3) 

where the subscript {6} refers to the removal of multiple bonds from a set of non-adjacent vertices {v} of 
the K p subgraph. Performing the HEC-typc homcomorphic expansion as above leads to the family 

HECk 1 ~2,k 2 -2,...,k p -2[(K p ){ b y.{ v } +G r ] (2.4) 



A special case of eq.(2.4) applies if in the starting family one removes all of the bonds connecting vertices 
in the K p subgraph to each other, so that this starting family is 

K P + G r (2.5) 

The H EC-type homeomorphic expansion of this family is thus 

HEC kl -2,k 2 -2,...,k p -2(Kp + G r ) (2.6) 

In this case one clearly obtains the same homeomorphically expanded family if one permutes the choices of 
kj: 

HECk 1 -2,k 2 -2,...,k p -2{K p + G r ) — i?£'C 7r ( fel )_2 

(K p + G r ) (2.7) 

where ir is an element of the permutation group onp objects, S p . The number of vertices in the homeomorphic 
classes of families (2.2), (2.4), and (2.6) is the same, namely 

v(HEC kl -2M-2,...,k p -2[{Kp) h + G r }) = v(HEC kl -2M-2,-,k p -2[(Kp) {b} . {v} + G r ] 



lHECk 1 -2,k 2 -2,...,k p -2[K p + G r ] 
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= p + r(l-2p + 5> J ) (2.8) 

3=1 

Note that for p = 2, since one can only remove b = 1 bond in the K 2 subgraph, thereby obtaining (^2)5=1 = 
K 2 , it follows also that 

HEC kl - 2M - 2 [(K 2 ) b=1 + G r ] = HEC kl -2M-2\K 2 + G r ] (2.9) 

For the case G r = K r , the chromatic number is 

x(HEC kl - 2M - 2 ,..., kp - 2 \Kp + Kr\) =2 (2.10) 
For G r =T r , for the cases r > 2 of interest here, 

X (HEC kl -2te-2,...,k P -2\K v + T r ]) = 3 (2.11) 

so that for r > 2, 

P(HEC kl - 2 ,k a -2,...,k p -2\K p +T r ],q = 2) = (2.12) 
For the r = 1 case, this family degenerates to a tree graph, 

p 

HEC kl - 2 ,k 2 -2,.,k P -2[Kp + Ti] = T Vt where u t = 1 + ]T(fc, - 1) (2.13) 

i=i 

with chromatic polynomial P(T Vt , g) = g(q — l)"* -1 and chromatic number \ = 2. 
A general result for the chromatic polynomial in the case r = 2, G r = T r is 

P 

p(ffEC fcl _2 lfc!l -2,...,fc p -2[X , J ,+T2],«) =«(«- 1)11^-1 (2-14) 



Since the respective indices 2kj — 1 of each of the D 2kj -i in eq. ( 2.14 ) are odd and since 

D kodd = (q-2)xpol(q) (2.15) 

where pol(q) is a polynomial in q of degree k — 3, it follows that 

p(if£C fel _2, fc2 -2,... )fep -2[# P + T 2 ],«) = <?(<? - l)(q - 2)f x poi'(g) (2.16) 

where pol'(q) is a polynomial of degree 2 X^=i(% — 2). 

Since the number of vertices n = v is a linear function of p, r, and fci, ...,k p , there are several ways of 
producing the limit n — + 00 (L denotes limit): 

L p : p — ► 00 with r ,k\ k p fixed (2-17) 

L r : r — > 00 with p,ki /c p fixed (2-18) 

ifcj : fcj — > 00 with p ,k\ , .., fcj_i, fcj+i, ...fc p fixed (2-19) 
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and, for the case where all of the kj's are the same, i.e., kj — fc c & Vj (where the subscript cb denotes 
"connecting bond"), 

Lk '■ k c t> — > oo with p 1 r fixed (2.20) 

As discussed before fl3| l , the limit L p is not of much interest from the viewpoint of either statistical mechanics 
or graph theory, since for any given graph G r and for any given (finite) value of q S Z+, as p becomes 
sufficiently large, one will not be able to color the graph (K p )b + G r or homeomorphic expansions thereof, 
and the chromatic polynomial will vanish. The limits L kj and Lk are also not of primary interest here since 
they generically yield compact boundaries B, as will be illustrated below. We shall therefore concentrate on 
the limit L r . The lowest choice of p that yields a noncompact boundary B(q) is p = 2. 

For general p and G r , we remark on two special types of HEC homeomorphic expansions. A minimal 
case is one in which the above homeomorphic expansion is carried out only on the bonds connecting a single 
vertex, taken to be vi with no loss of generality, in K p to vertices in G r , so k 2 = k 3 = ... = k p = 2: 

HEC kl -2,o„...,o(K p + G r ) (2.21) 

A particularly symmetric case is the one in which the homeomorphic expansions are the same on all of the 
connecting bonds: 

HEC kl -2M-2,...,k p -2{Kp + G r ) with fei = k 2 = = k p = k cb (2.22) 

In this case, the right-hand side of eq. ( fO| ) reduces to v = p + r +pr(k c b — 2). We proceed to give detailed 
analyses of some special families of homeomorphic expansions. 



III. CLASS OF FAMILIES HEC kl -i M - 2 (K 2 + K r ) 



The family HEC^-vfo-iiKz + K r ) has the special property that the graphs of this family only depend 
on the sum of ki and k 2 , not on each of these parameters individually: 

HECk 1 -2M-2(K2 + Kr)=HEC k > 1 -2,k' 2 -2(K2 + Kr) «■ h + k 2 = k[ + k' 2 (3.1) 

To incorporate this symmetry, we define 

Hk.r = HECk- 3 (K 2 + K r ) = HEC kl -2M-2{K2 + Kr) where k = ki + k 2 - 1 (3.2) 

As will be evident in our explicit results to be given below, the fact that this family has a noncompact 
W boundary B(q) in the L r limit can be understood to follow from its construction as a homeomorphic 
expansion of the starting family (^2)6=1 + G r = K2 + G r which (as was shown previously [Jk|) has a 
noncompact B (here, G r = K r ). The noncompactness of the boundary B obviously implies that in the L r 
limit, the chromatic zeros have unbounded magnitudes (in the q plane), since B arises via the coalescence of 
these chromatic zeros in this limit [Of R We note that the case k% = k2 = 2, i.e., k = 3 is a special case of 



4 For a family o f grap hs Gk, r , rather than analyzing the continuous locus B resulting from the L r limit in eq. ( 2.1E ) 
or the Lk in eq. ( |2.20[ ) , one may formulate a different problem: let {qo}k,r denote the set of chromatic zeros of Gk.r, 
and consider the union of this set, summed over both r and k, denoted {gojvfc.r = £/r=i X^fc^r^ojfc.r- This problem 
has been considered for the family Gk,r = Hk, r by A. Sokal, who finds (private communication) that the image of 
this union is dense in the vicinity of the origin, y = z = 0. 
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the graphs with noncompact B(q) that we have constructed and studied before |lq| , so we concentrate on 
the homeomorphic expansions k > 3 here. The number of vertices is given by the p = 2 special case of eq. 



(2.8) with ( |3.2| ), viz., v{Hk. r ) = (k — 2)r + 2. For arbitrary k > 2 and r > 1, -Hfc jr is bipartite, i.e., 

X(flfc.r) = 2 (3.3) 

For r > 2, the girth is g(H k:r ) = 2fc — 2. Indeed, all (non-self-intersecting) closed paths are of this length. 
For r = 1 and r = 2, the family degenerates according to 



and 



2k-2 



(3.4) 



(3.5) 



where, as defined above, T n and C' n are, respectively, the tree and circuit graphs with n vertices. Hence, we 
shall take r > 3. 

By use of the deletion-contraction theorem ^ we obtain the recursion relation 



P(H kir ,q) = D k P(H k , r _ u q) + {-l^qiiq - lJDfc-if- 1 
which we solve to get the chromatic polynomial 



P{H Kr ,q) = q{q-l) 



D 2k -2(D k y- 2 - (q- l)(D k ^ f (D k ) r - 2 -{(q- l)D k ^} r 



This has the form of eq. (L2j) with N a = 2, 



oi = Di 



and 



a 2 = (q- l)-Dfc-i 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



In the limit L r of eq. (2.18), the nonanalytic boundary locus B is determined as the solution of the degeneracy 
of magnitudes of the leading terms 



02 



(3.10) 



To determine B, we multiply eq. ( |3.1C| ) by |(?(t?-l)| to obtain |P(C fc ,g)| = l)P(C fc _ 1 , g)| (the resultant 
spurious solutions thereby introduced at q = 0,1, i.e., y = — l,oo, are understood to be ignored in the 
following discussion). Dividing by \a\ k yields 



5 We recall the statement of the addition-contraction theorem Q]- let G be a graph, and let v and v' be two 
non-adjacent vertices in G. Form (i) the graph G a dd. in which one adds a bond connecting v and v' , and (ii) the 
graph G con tr. in which one identifies v and v' . Then the chromatic polynomial for coloring G with q colors, P(G,q), 
satisfies P(G, q) — P(G a dd. , q) + P(G con tr.,q)- The reverse process of bond deletion, leading to the same equation, is 
known as the deletion-contraction theorem. 
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|i + (-i)V- 1 l = |i + (-i)*~V _a l (3-11) 

Since y = is a solution of eq. ( [3.11 ), B is noncompact in the q plane, passing through z — 1/q = 0, as noted 
above (equivalently, eq. (3.10) satisfies the condition given in the theorem of section IV of Ref. jl3| that 
guarantees that the solution locus B is noncompact in the q plane). (In contrast, [for fixed k] B is compact 
in the z or y plane.) 

Furthermore, we find that (i) 2{k — 2) curves on B, consisting of [k — 2) complex-conjugate (c.c.) pairs, 
intersect at z = y = 0; (ii) these curves approach the point z = y = at the angles 



(2j + 1)tt 



j = 0,...,fe-3 



(3.12) 



2(fc-2) ' 

Thus, if and only if k is odd, these angles include ±7r/2, i.e., a branch of B crosses the point z = y = 
vertically. To prove these results, we write (3.11) in terms of polar coordinates using eq. (1.5): 



Js-2 



p k - 2 (p 2 - 1) + 2(-l) fc { ( ocos[(fc - 1)0 + cos[(fc - 2)/3]} 



= 



(3.13) 



As p — + 0, so that z — > y (whence /3 — > 0), the solution is given by cos[(fc — 2)/3] = 0, which yields the results 
(i) and (ii) . In the terminology of algebraic geometry [ p0[ , the point z = y = is a singular (multiple) point 
on the algebraic curve B of index k — 2. 

As a further result, we find that (iii) if k is odd, B crosses the real q axis once, at a value q c that increases 
monotonically from 3/2 for k = 3, with limfe—^ g c = 2; (iv) if k is even, B never crosses the real q axis, so 
q c is not defined. To show these results, we first set (3 — in eq. (3.13), which becomes 



fc-3 



k-2 



s=0 



= 



(3.14) 



For odd fc, this has a single real, positive solution for p, which decreases monotonically from p = 2 for k = 3, 
approaching p = 1 as k — > oo. For even fc, eq. (3.14) has no (real) solution for p. For /3 = ir, eq. ( |3 .13 ) yields 



fc-3 



Jc-2 



2^> =0 



(3.15) 



Aside from the root p = 1, i.e., y = — 1, which is spurious as noted above, eq. ( 3.1 5| ) has no real roots. The 
result for q follows directly. 
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2.0 



1.0 



0.0 



-1.0 




(c) 



-2.0 



-2.0 



-1.0 0.0 1.0 

Re(z) 



2.0 



3.0 



FIG. 2. Boundary B in the z = 1/q plane for lim r ^oo H^^ with k = (a) 4 (b) 5 (c) 6. Chromatic zeros for H^^ r with 
(k,r) = (a) (4,30) (b) (5,18) (c) (6,10) are shown for comparison. 
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2.0 



1.0 



0.0 



-1.0 




-2.0 



-2.0 -1.0 0.0 1.0 

Re(y) 



2.0 



3.0 



FIG. 3. Boundary B in the y = X/(q — 1) plane for lim, — >cx) r with k = 6. Chromatic zeros for r with r = 10 are 
shown for comparison. 



The curves comprising B divide the z (equivalently, y or q) plane into k — 1 regions. For even k = 2£ 
and odd k = 21+1, these include (£ — 1) complex conjugate pairs of regions. In the q plane, B consists of 
k — 2 disjoint curves (a line for k = 3); each curve extends inward from complex infinity, turns around and 
heads back out to complex infinity along a different direction. For k = 3, B is a circle in both the y and z 
planes, given by the equation \y — 1| = 1, i.e., \z — h\ = 5- For k — 4, 5, 6 we show our calculations of B in 
the z planes in Fig. The corresponding plots in the y plane look more like (slightly distorted) cloverleafs, 



reflecting the simplicity of the degeneracy equation (3.11) in the y variable; as an example, we show the case 
k = 6 in Fig. ^. (Additional figures in the y plane are given in Ref. |2l]]). As is clear from Fig. |^, for k > 4, 
B has support for Re(q) < (equivalent to Re(z) < 0). Note that Re(z) = does not, in general, imply 
that z — 0, since for k > 4, B intersects the imaginary z axis away from the origin, e.g., at z = ±Zi, where 



Zi 



0.632456i for k = 4 



(3.16) 



Zi = iy/8/7 = 1.06904i for k = 5 



(3.17) 



(i/3)[10 T %/82] 1/2 = 0.32397H, 1.45508i for k = 6 



(3.18) 



In the y plane, the intercepts y = ±ipi of B with the imaginary y axis (aside from y = 0) are given, for even 
k = 21, by 
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pf- 1 \ P ^-l)-2(-lY = (3.19) 



and, for odd fc = 21 + 1, by eq. ( 3.1£ ) multiplied by p. Solving these equations for general k, we find that 
(with fc > 3) (i) for k — (0 or 1) mod 4, pi is nonzero, decreasing monotonically from \/2 for k = 4, 5 toward 
1 as fc — > oo; (ii) for fc = (2 or 3) mod 4, there is no nonzero intercept pi. Although B extends to Re(z) < 
(equivalently, Re(q) < 0), it never includes support for negative real z or q. 

In Fig. U we also show illustrative chromatic zeros for (a) k = 4, r = 30 (=> n = 62), (b) fc = 5, r = 18 
(=> n = 56), and (c) k = 6, r = 10 (=> n = 42). Aside from the always-present zeros at q = 0, 1 (z = oo, 1), 
the chromatic zeros generally lie close to the asymptotic curves comprising B onto which they coalesce as 
fc — * oo. For a given set of (finite) fc and r, the moduli of the zeros are bounded, and hence they avoid the 
points z = y = 0, approaching these only as r — > oo for a given k. 



As before [12 13], we let R\ denote the region including the positive real q axis for q > q c . For fc even, 
this includes the entire real q axis. For odd fc, we denote the region containing the rest of the real q axis as 
Ri,k odd] from our statement above, it follows that there is no Ri phase for fc even. For q € Ri, we calculate 

W([ lim H k . r ],q;qe R 1 ) = ( Y 2 ¥tvtl -l 1 - S ° dd ( 3 ' 2 °) 
v rwoo J |_ [ai(q)\ /[ ' ii k is even 

As discussed in Ref. |l^] , in regions other than Ri , there is no canonical choice of phase in taking the 1 / n'th 
root in eq. (|j~l|), so that one can only determine |VT({G}, q)\. We find 

|W([lim H k>r ],q)\ = Mg)| 1/(fc-2) , for fc odd and q £ R 2 . k odd (3.21) 

r— >oo 

|T^([lim i/ fc , r ],g)| - |a 2 (Q')| 1/(fc " 2) , for fc odd and q # R u R 2 . k odd 

and for fc even and q $ i?i (3.22) 

On B, \W\ is continuous but nonanalytic. 



IV. CLASS OF FAMILIES HEC kl - 2 ,o,o(K 3 + K r 



When p > 3 in the homeomorphic class of families ( [2.2] ) or (2.6) with G r = K r , the graphs and their 
chromatic polynomials depend on the individual values of the kj, not just on the sum, as in eq. ( |3.l| ). A 



simple homeomorphic class of families of this type is the special case of eq. ( 2.21 ) for p = 3 and G r = K r , 
namely 

6 fc , r =HEC kl - 2i o, (K 3 + K r ) , h = k (4.1) 
We have u(Ofe, r ) = r(k — 1) + 3 and 

X(e fc , r ) = 2 (4.2) 
For our present purposes, it will suffice to study the family 03, r . We calculate the chromatic polynomial 
P(&3,r, q) = (q- 2) 2j - 4 P(e 3 , 2 , g) - q(q - 1)(D 4 ) 2 \(q - 2) 2 ^ 4 - (D^ 
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-q(q - l)(g - 2)(q 2 - 4q + 5) 2 [(q - 2) 2r ^ - (q 2 - Aq + 5)^ 2 

where 

P(®3,2, q) = q(q - 1) [(« - l)Ai - (« ~ 2)D 6 
This has the form of eq. ( |1.2| ) with iV = 4 and 

ai = -D4 = g 2 — 3g + 3 

a-2 = q 2 - 4<7 + 5 



and 



03 = (?-l)(g-2) 



a 4 = (q- 2) 2 



(4.3) 

(4.4) 

(4.5) 
(4.6) 
(4.7) 

(4.8) 



0.6 




-0.6 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-0.5 0.0 0.5 1.0 

Re(z) 

FIG. 4. Boundary B in the z = l/q plane for lim T ._ >00 03, r . Chromatic zeros for 03, r with r = 16 are shown for comparison. 
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The nonanalytic boundary B for the limit as r — ► oo is shown in Fig. |4|. As with the Hh.r family discussed 
above, and the other families to be discussed below, one sees the important feature that B includes support 
for Re(z) < or equivalently Re(q) < 0. The boundary divides the z (or equivalently y or q) plane into six 
different regions: (i) i?i, including the interval < z < 1/2 on the positive real z axis; (ii) R2, including 
the union of the intervals z > 1/2 and z < 0; (iii) a complex-conjugate pair of regions R% and R% lying just 
above and below R\ in the z plane; and (iv) a second pair of complex-conjugate regions R± and R\ located 
such that i? 4 lies above and adjacent to R 3 (hence between R 3 and R 2 ) and R± lies below R^. In the regions 
Ri and R%, a\ and a,2 are dominant, respectively, while aj is dominant in regions Rj and R* for j ' = 3 and 
j = 4, respectively. We have 

W([lim e kl r],q) = (ai) 1/2 for q e Ri (4.9) 

r — »oo 



\W([ lim fc , r ],g)| = \a 2 \ 1/2 for g € i? 2 (4.10) 

r — >oo 

and 

\W([lim e k , r ],q)\ =\aj\ 1/2 for q G R h m ,j = 3, 4 (4.11) 

r— >oo 

The region boundaries between regions i?i and i?j are the solutions of the respective degeneracy equation 
| Oil = I where a,i and aj are leading terms. Dividing these equations by \q\ 2 and re-expressing them in terms 
of z to get the corresponding degeneracy equations in the z plane, one has |5j | = \a,j | with a\ = 1 — 3z + 3z 2 , 
a,2 = 1 — 4z + 5z 2 , 5a = (1 — z)(l — 2z), and 04 = (1 — 2z) 2 . Clearly, each pair of the degeneracy 
equations |aj| = \dbj\ in the z plane has a solution at z = 0, which shows that, in accordance with the 
condition of Ref. p3[ , S is noncompact in the g plane, passing through the point z = 0. Thus, as is evident 
in Fig. ||, all of the regions are contiguous at z = 0, where six curves (i.e., three branches) of B meet in a 
multiple point. For small £ = |z|, the degeneracy equation for the part of B separating R\ and R 3 has the 
form 2£ 2 (1 — 2 cos 2 9) + 0(C 3 ) = 0, so that this boundary crosses the point z = at the angles given by 
cos0 = 1/V2, i.e., 

9 Ri ,r 3 , 0r uR *=±^ (4.12) 

(The other solution, cos6> = —l/y/2 is not relevant because in this region, a\ and az are not leading terms.) 
Similarly, for small £, the degeneracy equation for the part of B separating the regions R3 and R4 is 
£(3£ — 2 cos 6) = so that as C — > 0, this boundary crosses the point z — at the angles given by cos 9 = 0, 
i.e., 

0R 3 , Ri , 0fl 3)fl j=±| (4.13) 

The same type of reasoning applied to the i?2, R4 and i?2, -R4 boundaries shows that they cross z = at the 
angles 

37T 

9a,« 41 9a ) flj=± T (4-14) 
Concerning the boundary separating regions R\ and R2, we observe that the relevant degeneracy equation 

B{R 1 ~R 2 ): |l-3z + 3z 2 | = |l-4z + 5z 2 | (4.15) 
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has the solution z = 1/2 as its only solution other than z — where eii and a-i are leading terms. Hence, 
q c = 2 for this family. Concerning the two complex-conjugate (c.c.) pairs of multiple points, we note that 
the multiple point 22,3,4 where regions R2, R3, and i?4 are contiguous, and its c.c, are 

*2,3,4, < 3 ,4 = S-V^i^/e = 1 ± _i_ (4 . 16) 

Hence this c.c. pair lies on the unit circle \y \ = 1 or equivalently \q — 1| = 1 in the respective y and q planes: 

2/2,3,4, 2^3,4 = e ±W3 (4.17) 

92,3,4, 9 2 *,3,4 = i 1/2 e^l & (4.18) 

A corresponding analysis can be given for the multiple point zi.2,3 where regions R\, R2, and R3 are 
contiguous. 



V. CLASS OF FAMILIES HEC kl -2M-2(K 2 + T r ) 



Illustrative graphs of the family 

HEC kl - 2M -2(K2+T r ) (5.1) 

with (fei, k 2 ) = (3, 2) and (fci, £2) = (3, 3) are shown in Fig. |](c,d), respectively. From the general equation 
(2.8) with G r = T r and p = 2, we have 

v(HEC kl - 2M - 2 (K2 + T r )) = 2 + r(fci + fe 2 - 3) (5.2) 

The special case k\ = = 2 is the family i^2 + T! r = (i^2)fc=i + 2V which we analyzed in detail earlier so 
we concentrate on the homeomorphic expansions of this starting family here, i.e., k\ > 3 and/or &2 > 3. In 
order to calculate the chromatic polynomial for the graphs in this homeomorphic class of families, we shall 
utilize a generating function method as we did in Ref. Jl(f | for infinitely long, finite-width strip graphs of 
various lattices. The generating function T is a function of a symbolic variable x and yields the chromatic 
polynomials P[HECk 1 -2,k 2 -2(K 2 + T r ), q ) via the Taylor series expansion around x = according to 



00 

T^HEC kl -2M-2(K2+T r ),q,x^ =J2 P { HECk i- 2 > k ^(K2 + T r ),q^)x r (5.3) 

r=l 

The summation starts at r — 1 since this is the minimum value of r in this family. As in our earlier work 
Jlql, we find that T is a rational function of x and, separately, of q, of the form 



, _ \ Af[HEC kl -2,k 2 -2(K2+T r ) iq ,x) 

T{HEC kl -2M-2{K2 + T r ),q,x) =— } f (5.4) 

V ' v(HEC kl -2, k2 -2(K2 + T r ),q,^ 



with 

N(HEC kl -2M-2(K2 + T r ),q,x) =xJ2 A M xj ( 5 - 5 ) 
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(so that deg x {J\f) = jmax + 1; the prefactor of x reflects the fact that the minimum value of r is 1) and 



Jmax. 

v(HEC kl - 2M - 2 (K 2 + T r ),q,x^ = 1+ >T b r {q)x j ' 



(5.6) 



where j max and j' max depend on the specific family, the Aj and by are polynomials in q, and their dependence 
on k\ and k 2 is left implicit. (This notation follows that of Ref. |l6j except that in eq. (pli|), we use A rather 
than a to avoid confusion with the aj functions in eq. fll.2p .) A general formula for A$ is 

A = P(T kl+k2 _ u q) = q(q - (5.7) 

The other Aj's and bj> polynomials will be given below for specific families. 

For the starting family, HECq$(K 2 + T r ) = K 2 + T r = (K 2 )b = i + T r , and indeed its generalization 
(K p ) b + T r , one has @ 



P((K p ) b + T r , q) = q^ Uq-p- l)^ 1 + b(q - p) 



(5.8) 



(recall the notation in eq. (1.8)). The equivalent representation in terms of a generating function is obtained 

(5.9) 



by noting that eq. (5.8) has the form of eq. (1.2) with cq = and N a — 2, i.e., 

P((K p ) b + T r ,q) = Cl ( ai ) r + c 2 {a 2 y 
with 

ai = q - p - 1 , a 2 = q -p 
„(p+i) 



ci 



and 



c 2 = 



q-p-1 

bq (P+i) 



(5.10) 
(5.11) 

(5.12) 



{q - p) 2 

We find that the chromatic polynomial P((K p )b + T r> q) is given as the coefficient of x r in the Taylor series 
expansion about x — of the generating function 

x(A + A x x) 



r 



(<V \ , rr \ X[A +A 1 X 



where 



Xj = aj , j = 1,2 



A Q = aiC\ + a 2 c 2 



and 



Ai 



-aia 2 (a + c 2 ) 



(5.13) 



(5.14) 



(5.15) 



(5.16) 



For a general graph G, this type of relation between the form for the chromatic polynomial (5.9) and the 



generating function (5.13) has a straightforward generalization to the case where N a > 3. In the present 
case (-^2)5=1 + T r = HECo o(K 2 + T r ), we find, for the upper limits on the sums in Af(x) and T>(x) the 
values j max = 1 and j' max = 2. 
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A. Families HEC kl - 2 ,o(K 2 + T r ) 



We first consider the simplest homeomorphic expansion of type HEC on the starting family K 2 + T r , in 
which one performs this expansion only on the bonds connecting one of the two vertices in the K 2 subgraph 
with the vertices of the T r subgraph. With no loss of generality, we choose this vertex in K 2 to be v%, so 
that k 2 — 2 = and k x — 2 > 1. For brevity of notation, define 



Kr = HEC kl ^ 2fi (K 2 +T r ), k x = k 



(5.17) 



with v(Ok, r ) = 2 + r{k — 1). When referring to the collection of the graphs of this family for various r, we 
denote it as {Ok}- The case k = 2 is just the original starting family, K 2 + T r . For the family Ok, r with 
k > 3 we find that 



jmax — 2 , imax ~ ^ 

For general k > 3, the denominator of the generating function can be written as 

V{{O k }, q, x)=\l-(q- 2)D k x] \l - (q - 3)D k x - (q - l)(q - 2)D k _ 1 D k x 2 



(5.18) 



where 



and 



(1 - X 2 x)(l - Ai p a;)(l - Xi m x) 



(q - 3)D k ± {(q - 3)D k } 2 + 4{q - l)(q - 2)D k _ x D k 



1/2 



(5.19) 



(5.20) 



\ 2 = (q-2)D k (5.21) 
The coefficient functions in Af are A , given by the special case of eq. (|5.7|), viz., A = q(q — l) k , and 



M = -q(q-l)D k 



(2q - 5)(q - I)*" 1 -(q-2) qD k - 2(-l 



(5.22) 



A 2 = -q 2 {q-l) 2 {q-2)D k ^Dl (5.23) 
(In deriving the expression for A X) we have used the identity D 2k ^i = [qD k — 2(— \) k ]D k .) We calculate 

W({O k },q) — (\max) 1 ^ k ~ 1 ' ) for qeR, (5.24) 

and 

\W({O k },q)\ = \X max \ 1 ^ k - 1 '> for q e R, + R x (5.25) 



where A maa; (g) denotes the leading A in eq. (5.19) in the respective regions of the q plane. For a general 
family {G} of graphs with a generating function 



T({G},q,x) = ^ = P(G r ,q)x r 

3=30 



(5.26) 
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with 



= n a - a, 

j'=l 



(5.27) 



the boundary B in the limit r — > oo is given (see eq. (4.16) of Ref. |na) as solution locus of the equa- 



tion expressing the degeneracy in magnitude of two leading A terms in eq. ( 5.27 ), where W switches form 
nonanalytically, 



|A mQX (g)| = |A 

ma x- (q)\ 



(5.28) 



As an illustration, we consider the lowest homeomorphic expansion, k = 3, for which eqs. ( 5.2C| ) and (5.21) 
read 



A 



lp,m — 2 V? 2) 



(g - 3) ± (g 2 - 2q + 5 



1/2' 



(5.29) 



A 2 = (g - 2) 2 



(5.30) 



Let us now write Ai p , Ai m and A2 in terms of the variable z = 1/q in polar coordinates, as given by eq. ( |1.4| ) 
and Taylor-expand the resulting expressions for small £. We obtain for the magnitudes squared of the A's 



|Ai p | 2 = 1 -4Ccos6> + 2C 2 (cos(26>) + 2) + 0(C 3 ) 



(5.31) 



|A lm | 2 = c 2 + o(C 3 ) 



(5.32) 



|A 2 | 2 = l-4Ccos0 + 4C 2 + O(C 3 ) 



(5.33) 



Thus, for small values of \z\, the boundary B is given by the equation | Ai p | = | A2 1 , which yields ( 2 cos(20) + 
0(C 3 ) = as C — ► 0. Hence, B is noncompact in the q plane and crosses the origin of the z plane at angles 
such that cos(2#) = 0, i.e., 



7T 7T 

— +n— 
4 2 



< n < 3 



(5.34) 



This is evident in Fig. [|. From eq. (1.3), it follows that y — > z as z — > 0, so that these angles are the same in 
the y plane. The boundary B divides the z plane into four regions: (i) R±, including the interval < z < 1/3 
on the real z axis; (ii) R2, including the intervals — 00 < z < and 1/3 < z < 00 on the real z axis; and (iii) 
the complex-conjugate pair of regions R3 and R% lying roughly above and below z = 0. 
We find 



W([lim 3l r],g) = (Ai P ) 1/2 for q G ifc 



(5.35) 



|W([ lim Q3, r ],g)| = lAipl 1 / 2 for g e i?, 2 



(5.36) 



and 



18 



\W([\im 3 , r ],q)\ = \\2\ 1/2 for q e R 3 , R* 3 



(5.37) 



A portion of B crosses the positive real axis at z = z c = 1/3 so that q c = 3 for lim,,—^ 3r . Along 
this portion of the boundary |Ai m | becomes degenerate with |Ai p |, although the former never dominates in 
magnitude over the latter. This portion of the boundary ends in two c.c. multiple points. Note that 



W([lim 3r ],q) = at q = 2 

r— >oc 

consistent with eqs. ( Fill) and ( jyg) 



(5.38) 



0.6 



0.3 



0.0 



-0.3 



-0.6 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-0.5 0.0 0.5 1.0 

Re(z) 

FIG. 5. Boundary B in the z = 1/q plane for lim, — too Ok r with k = 3. Chromatic zeros for 0^ r (k, r) = (3, 29) are shown 
for comparison. 



B. Families HEC kl - 2 ,k 2 -2(K 2 +T r ) with k\ = k 2 

We next consider the case of symmetric homeomorphic expansion of the bonds from the two vertices of 
the K2 subgraph, i.e., k\ = = k. For brevity of notation, we denote 

U k>r = HEC kl _ 2M _ 2 (K 2 + T r ) , fei - k 2 = k (5.39) 

with v(Uk, r ) — 2 + r(2k — 3). When referring to the collection of the graphs of this family for various r, we 
denote it as {[/&}. 

For the family Uk.r we find that 
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jmax — 3 , jmax ^ 

The denominator of the generating function can be written as 



(5.40) 



V({U k }, q,x)= l-(q- 2)Dlx - (q - l) 3 D 2 



1 - D k (D k+1 - D k )x -{q- IfD^xDlx 2 



(5.41) 



Observe that each factor of x in T>({U k }, q, x) occurs with at least one accompanying factor of D k . We have 



^({Ukj^q, x) = (1 - \i p x)(l - Ai m a;)(l - A 2p a;)(l - A 2m x) 



where 



D,, 



> lp,m 



(q - 2)D k ± [{q - 2)D k } 2 + 4(q - l) 3 (Z? fe _x) ; 



1/2 



and 



*2p,m 



(D k+1 - D k ) ± {D k+1 - D k f + A(q - l) 2 DnA 



1/2 



For the polynomials in the numerator Af(x), eq. (pj) gives Aq = q(q — l) 2k 2 and we find 



A 1 =q(q-1) 



(D 2k -i) - (I - 1) (? - 3)^fc + D k+1 D 



A 2 = -q(q - l) 2 D k _ x Dl [q 2 D 2 - q(q + 2){-l) k D k + 2q-\ 
q(q - 2){q - \) k ~ 2 D k + (q - 2){-\) k (q - \) k ~ 2 + 1 



We calculate 



and 



A^-q^q-lfDUDl 



W{{U k },q) = [\ max {q)] l K 2k -V for q G R t 



\W({U k },q)\ = \\ max (q)\ 1 /( 2k -V for q G Rj + R x 



(5.42) 



(5.43) 



(5.44) 



(5.45) 



(5.46) 



(5.47) 



(5.48) 



(5.49) 



where X max (q) refers to the respective dominant A in the given region. The boundary B is given, as before, 
by the solution locus of eq. ( 5.28| ). We show this boundary in Fig. ||(a), together with illustrative chromatic 
zeros for the case r = 19. From the property ( ^.15| ), it follows that the A's are degenerate at zero when q = 2, 



Aip, m (<? = 2) = \2p,m(q = 2 ) = for k odd 
so that B passes through the point q — 2 if k is odd, and, furthermore, 



(5.50) 
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W([ lim U k , r ],q = 2) = for k odd 



(5.51) 



This is consistent with the facts that x({^fe}) = 3 and P(Uk, r , q = 2) = as special cases of eqs. ( 2.11 ) and 



(M3- 

We next give some explicit results for the lowest two cases, k = 3 and k = 4. For the family Uk=3, r eqs. 



(5.43) and (5.44) yield 



Alp,m;fc=3 — TJ^ _ ^) 



(g - 2) 2 ± (V - Aq 3 + 12q 2 - 20q 
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1/2' 



(5.52) 



A2p,m;fc=3 — -^(q — 2) 



(<f - 4g + 5) ± g 4 - 4g J + lfV - 20g + 17 



1/2" 



(5.53) 



To investigate the boundary in the vicinity of z = 0, we divide the degeneracy equations by |g| 3 and express 



the results in terms of z = 1/q in polar coordinates, as given by eq. (1.4). We then Taylor-expand these 
equations for small £ = \z\. This yields 



|Alp;fc= 



4 - 24Ccos<9 + 4C 2 [8cos(260 + 9] - 8£ 3 [12cos6> + cos(30)] + 0{( 4 



(5.54) 



|A 2 p ; fe =3 | 2 = 4- 24Ccos6» + 4C 2 [8cos(26») + 9] - 16C 3 [6cos6> + cos(36>)] + 0(C 4 ) 



(5.55) 



|Ai m;fc=3 | 2 = 4C 2 + 0(C 3 



(5.56) 



|A 2m ; fc = 3 | 2 = 4C 2 + 0(C 3 ) 



(5.57) 



Thus, in the vicinity of z — the boundary B is given by the equation |Ai p; fe = 3| = | A2p ; fe=3 1 , which yields 
cos(30) = for £ ^ 0. Hence, six curves on B (forming three branches) cross the point z — (and hence also 
the point y = 0), at angles 



6 +n 3 ' 



< n < 5 



(5.58) 



The boundary B divides the z plane into six regions. As one traverses a circle around the origin, z = 0, 
starting with a small positive real z value, first Ai p; fe = 3 is dominant, and then the dominant A's alternate 



between Ai p; fc = 3 and A2 P ;fc=3- The resultant W functions are given by eqs. (|5.48| ) and (|5.49| ) 



For the family Uk=i, r eqs. (5.43) and (5.44) yield 



Aip, m; fc =4 = 7^{q - 2)(g 2 - 3g + 3) 



(q 2 - 3g + 3) ± (g 4 - 2q 3 + 3q 2 - 6q 



1/2 



(5.59) 



^2p,m;k=4 = -^(l — 3g + 3) 



(q 3 - 5q 2 + 9q- 7)± 



6g 5 + 15g 4 - 24g 3 + 35g 2 - 42g + 25 



1/2" 



As before, the boundary B is determined by eq. (5.28). This boundary is plotted in Fig. ra(b). Let us 



(5.60) 



write Aip, m and \2 P , m in terms of the variable z = 1/q in polar coordinates, as given by eq. (1.4) and 
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Taylor expand the resulting expressions for small £. We find that the two A's that are dominant near z = 0, 
namely \\ p and A2 P , have squared magnitudes that coincide through 0(C 4 ) and differ in order 0(C 5 ) in 
the coefficient of £ 5 003(56*). (The other two A's are subdominant near z — 0; indeed, they vanish there, 
as |Ai m | 2 ~ |A2m| 2 = 4£ 2 + 0(( 3 ).) It follows that as C — * 0, the boundary B is given by the equation 
|Ai p | = |A 2p |, which yields cos(5#) = for ( ^ 0. Hence, ten curves, forming five branches of B, cross the 
point z = (and also y = 0) at the angles 

= — +n~ , 0<n<9 (5.61) 
10 5 

Here, the boundary B divides the z plane into eight regions. As one traverses a circle around the origin 
starting with a small positive real value of z, first Ai p; fc = 4 is dominant, and then the dominant A's alternate 
between Ai p; fc = 4 and \2 P -,k=4, as before with k = 3. The resultant W functions again follow from eqs. (5.4S) 
and (5.49). Thus, comparing our results for Uk, r with k — 3 and fc = 4, we observe that the number of 
curves on B passing through z — 0, and the number of regions, increase as k is increased. One should also 
remark that certain properties of the region depend on whether k is even or odd, such as the fact that B 
passes through z = 1/2 for odd k > 3. 
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1.2 




-1.2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

Re(z) 

FIG. 6. Boundary B in the z = l/q plane for limr—.oo T with k = (a) 3 (b) 4. Chromatic zeros for r with (k,r) = (a) 
(3,19) (b) (4,15) are shown for comparison. 

VI. DISCUSSION 

In this section we discuss general features of the Potts antiferromagnet W functions on all of the families 
of graphs that we have constructed and studied in this paper and our earlier work, with loci B that are 
noncompact in the q plane. In (the theorem of section IV of) Ref. |l3| we gave the general algebraic 
condition that for a particular family of graphs, the infinite-vertex limit yields a locus B that passes through 
z = 0. From our calculations we have observed a geometrical regularity in the families of graphs that have 
this property, viz., that in this infinite- vertex limit they all contain an infinite number of different, non- 
overlapping (and non-self-intersecting) circuits, each of which passes through at least two fixed, nonadjacent 
vertices. This immediately implies that these aforementioned nonadjacent vertices have degrees A that go 
to infinity in this limit. One thus sees at the graphical level why the derivation of the large-g series for 
the reduced function W re d({G},q) or equivalently W({G},q) on regular lattices with free or some type of 
periodic boundary conditions works; in the thermodynamic limit, such lattices do not have the two or more 
nonadjacent vertices with degrees A — > oo as in the criterion stated above for noncompact B. 

We remark that in studies of the Potts model on regular lattices, it has been useful to utilize certain 
boundary conditions that do involve vertices which, in the thermodynamic limit, have infinite degree A, 
since these make it possible to preserve exact duality on finite lattices p^-|25|. On the square lattice, the 
duality-preserving boundary conditions (DBC's) of types DBC-1 and DBC-2, in the notation of Ref. |2~5}] 
involve one such vertex with A — > oo while the DBC-3 type involves two such vertices; however, in the latter 
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case, these are adjacent. Hence, none of these duality-preserving boundary conditions invalidates a large-g 
expansion of W re d{{G}, q) for this lattice. 

An elementary topological property should be noted: for an arbitrary family of graphs, the continuous 
locus of points B where W is nonanalytic in the infinite-vertex limit does not necessarily enclose regions in 
the q plane. Indeed, in Ref. Jl6[ we carried out exact calculations of a variety of infinitely long, finite-width 
strips of different lattices and found loci B that consisted of arcs (and line segments) that did not enclose 
regions. (We have also performed similar calculations for other types of strips with specific types of end- 
graphs that yield loci B that do enclose regions.) None of these loci were noncompact in the q plane. In 
contrast, if the infinite- vertex limit of a family of graphs yields a locus B that passes through z = 0, then, 
given that it is not a semi-infinite line segment on the positive or negative z plane, it necessarily separates 
the z and equivalently the q or y planes into at least two different regions where W is an analytic function. 
In Ref. Jl6| , we observed that the endpoints of the arcs (or line segments on the real q axis) that comprised 
the nonanalytic locus B were determined by the branch points of certain algebraic expressions occurring 
in the A's that entered in degeneracy equations of the form of eq. flS.28 ). Since, for example, the A's for 



the family {Uk}, eqs. (5.43) and (5.44), also contain square roots, one may wonder what role the branch 
points of these roots play with regard to the boundary B. We have investigated this and have found that 
the analogues of the arcs extending between branch points of the square roots, which comprised the various 
boundaries in Ref. [jl6|| , do not yield endpoint singularities on B here. The reason for this is that either (i) 
these arcs involve degeneracy of A's that are nonleading in a given region, or (ii) although they coincide on 
part of their length with the actual boundary, the portion of the arcs containing the endpoints lies off this 
boundary, in a region where condition (i) holds. This is illustrated in Fig. ^ for the {Uk} family with k = 3, 
where we show the actual boundary £?, as in Fig. |6](a), together with the arcs (drawn in a thicker black) 
formed by the degeneracy conditions 

|A 2p | = |A 2m | (6.1) 

and 

|Ai p | = |A lm | (6.2) 
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FIG. 7. Boundary B in the z = l/q plane for lim, — 100 Us,r, shown together with the locus of solutions of the degeneracy 



equations (J6 . l[) and (3.2). See text for discussion. 



The complex-conjugate pair of arcs which constitute the solution locus to eq. ( p.l[ ), with endpoints at 
z 2a, z^ a = 0.04106 ± 0.43627i and Z2b, z\ h = 0.54717 ± 0.08333i, lie in the interior of the regions R4 and 
i?4- But in these regions, neither A2 P nor \2 m is dominant; rather, as we discussed in the section on the 
Uk >r family, it is \\ p that is dominant in these regions. Hence, this complex-conjugate pair of arcs is not 
relevant to the actual boundary. As shown in Fig. [7], the locus of solutions for the degeneracy equation 



(5.2) forms a complex-conjugate pair of arcs that cross each other and the real z axis at z = 1/2 and 
have endpoints at z la , z\ a = 0.10169 ± 0.37529i and z u , z\ h = 0.73164 ± 0.12612i. The portion of these 
arcs that lie to the right of Re(z) — 1/2 are not relevant for determining the boundary because in this 
region, denoted R2 in our section above on the {Uk} family, the dominant A is A2 P . The portion of the 
complex-conjugate arcs lying between the multiple (crossing) point z — 1/2 on the right and the multiple 
points forming T-intersections on the left at Zt, z'^ — 0.3204 ± 0.2110z does coincide with the boundary, 
since on this portion the boundary is determined by the degeneracy of magnitudes of leading eigenvalues 
Aip and Ai m , i.e., eq. ( |6.2[ ). However, to the left of these T-intersection points, the arcs forming the solution 
locus of eq. (|6.2|) again deviate from the actual boundary, which is determined by the degeneracy of leading 
magnitudes |A2 P | = |Ai p |. Hence, once again, the left endpoints of these arcs are not relevant for B. Thus, 
the boundaries B for the families considered here, even aside from the fact that they automatically enclose 
regions owing to their noncompactness, in contrast to the arcs comprising the loci B for the strip graphs 
that we studied in Ref. [HJ, also differ qualitatively in that they do not contain any arc endpoints. We have 
explained why this is so even though some of the families do have A's containing branch points singularities. 
This emphasizes that the existence of the arc-like structure of the loci B that we found in Ref. |ll| depended 
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not just on the presence of branch point singularities in the A's entering via equations analogous to eq. ( 5.42|) 
in the denominators V of the generating functions for various strip graphs, but also on the fact that these 
branch point singularities occurred at endpoints of arcs in regions where these arcs were the solution loci of 
magnitude degeneracy equations for leading A's. 

From our studies of many homeomorphic expansions, we have found the general feature that for a given 
homeomorphic class parametrized by some set of homeomorphic expansion indices {kj}, the number of 
regions separated by the locus B in the r — > oo limit is a nondecreasing function of the above expansion 
indices. This is in accord with the constraints from algebraic geometry, in particular, the Harnack theorem 
pCfl . The application of this theorem is simplest in the case where B is a nonsingular algebraic curve (so that 
the number of regions N reg , = N comp . + 1 where N comp , denotes the number of connected components of B) , 
as is the case for the r — > oo limit of the families Tk <r — HEGt-2{K2 + T r ) and Su,r — HEGk-iiKa + T r ). 
In these cases (where B has no multiple points) the Harnack theorem states that the number of regions is 
bounded above by g + 1, where g is the genus of the algebraic equation in the variables Re{q) and Im(q) 
whose solution set is B, viz., g = (d— l)(d— 2)/2 where d is the homogeneous degree of this equation. Since, 
as our exact solutions show, this degree and the resultant genus increase as a function of the homeomorphic 
expansion indices, the Harnack upper bound also increases. It should be noted, however, that the number 
of regions may remain the same as one performs a homeomorphic expansion (e.g., N reg , — 2 for the r — > oo 
limit of HEG k - 2 {K 2 + T r ) for k = 2, 3, and 4, then N reg . = 3 for k = 5, 6, and 7, etc.) 

For algebraic curves B with singularities such as multiple points, one no longer has the simple relation 
N re g. = N comp . + 1. Indeed, in the families constructed and analyzed in this paper, B consists of a single 
connected component, i.e., N comp . — 1, while the number of regions N reg _ and the index of the multiple point 
at z — increase as a function of the indices {kj} describing the homeomorphic expansions. The origin of 
this increase is the increase in the degree of the polynomials in z that occur in the degeneracy equations 
relevant for the boundary B in the neighborhood of z = (e.g., eqs. ( 3.11 ) and the resultant ( 3.12 ) for Hk t r, 
etc.). 

The boundaries B found for the families of graphs in the present paper exhibit some interesting differences 
with respect to the (similarly noncompact) boundaries that we have previously studied |i~3( ] ; the latter shared 
three general properties: (i) a single curve on B passes through z — 0, so that this point is a regular point 
on B as an algebraic curve; (ii) if q € B, then Re(q) > 0, i.e. B includes support only for Re(q) > 0, or 
equivalently, for Re(z) > 0; and (hi) if z € B and Re(z) — 0, then z = 0. In contrast, for the families studied 
in the present paper, the homeomorphic expansions, which are all of the type HEC, yield, in the r — > oo 
limit, respective boundaries B that differ in each of these aspects: (i) several branches of the algebraic curve 
B pass through the point z — 0, which is thus a singular (multiple) point on this curve; (ii) B includes 
support for Re(q) < or equivalently Re{z) < 0, and (hi) B includes points with Re(z) — other than 
z = itself. Some insight into this can be gained by recalling the differences in the types of homeomorphic 
expansions of starting graphs. Since in the original starting families such as (K p )b + G r and K p + G r , it was 
the connecting bonds linking the vertices of the K p subgraph with the vertices of the G r subgraph that gave 
rise to the noncompactness of the respective boundaries B, it makes sense that homeomorphic expansions of 
these bonds would alter the nature of the boundaries at z — 0. Since this change involves multiple branches 
of B passing through this point at various angles, as in eqs. (frf), (fL34j ), (HH) and ( pT6l"l ), this also shows 



why this type of homeomorphic expansion leads to boundaries that include support for Re(q) < 0, but 
q ^ 0, oo, i.e., Re(z) < but z ^ oo, 0. In contrast, since the bonds within the G r and K p subgraphs, 
by themselves, are not directly responsible for the noncompactness of B, it is plausible that homeomorphic 
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expansions of these bonds would not change the nature of B at z = 0. 



VII. BOUNDARY FOR L k LIMIT 



Finally, we briefly discuss the boundary B that results when one takes the limit of eq. (2.20), i.e., 
k — > oo with r and p fixed, for the families studied in this paper. We consider the nontrivial range r > 2 
since for r — 1, families typically degenerate into tree graphs with B = 0. For the families studied here for 
which we have obtained the chromatic polynomials for general k, including Hk, r , Ok, r , and Uk,r, we find that 
this boundary is simply the unit circle \q — 1| = 1. This is easily understood, since one can reexpress our 



formulas for chromatic polynomials in terms of polynomials in a = q— 1 using eqs. (1.6) and (1.7). From eq 



(1.7), it follows that the chromatic polynomials are of the form of eq. (1.2) with a as the quantity entering 
in terms raised to respective powers proportional to k, together with possible Co terms involving (— 1) , and 
hence, as discussed in our earlier work fl2|| , B is determined by the condition \a\ = 1. Hence, q c = 2 and B 
has no support for Re{q) < 0. Since the locus \q — 1| = 1 is compact in the q plane and since our focus here 
is on the situation where B(q) is noncompact in the q plane, passing through z = 0, and the connection with 
the existence of large-g series for W re d., the limit is not of primary interest here. We mention, however, 
that B divides the q plane into two regions, {Ri)i, k and [B,2)L k , which are the exterior and interior of the 
unit circle \q — 1| = 1 (the subscripts are appended to distinguish these regions from those for the L r 
limit). It is straightforward to calculate the W function in these two regions in the Lk limit of the various 
families. For example, for the Hk, r family (for fixed r > 2), 



and 



W([\im H H , T ],q) = q-l for q G (Ri)L k (7.1) 

k— *oo 



\W([ lim H k , r ],q)\ = 1 for q € (R 2 ) Lk (7.2) 

k—too 



VIII. CONCLUSIONS 

In this paper we have presented exact calculations of the zero-temperature partition function, Z(G, q,T = 
0), and ground-state degeneracy (per site), W({G},q), for the g-state Potts antiferromagnet on a number 
of families of graphs {G} for which the boundary B of regions of analyticity of W in the complex q plane 
is noncompact and has the properties that (i) in the z = 1/q plane, the point z — is a multiple point on 
B; (ii) B includes support for Re(z) < 0; and (iii) B crosses the imaginary z axis away from z = 0. Our 
results yield insight into the conditions which preclude the existence of large-g Taylor series expansions for 
the reduced function W re d. — q _1 W . This insight is valuable since large-g expansions, where they exist, are 
of great utility in obtaining approximate values of the exponent of the ground state entropy, W. 
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